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ABSTRACT 

We present sensitive sub-mm imaging of the Seyfert 1 galaxy NGC 7469 
at 850 //m and 450 /xm with the SubmiUimetre Common User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) and ^^CO 
J=3-2 hne observations of its central starbursting region. The global dust 
spectrum, as constrained by the new set of sub-mm data and available 1.30 mm 
and IRAS 100 /im, 60 /im data reveals a dominant warm dust component 
with a temperature of Tj ~ 35 K and a global molecular gas-to-dust ratio 
M(H2)/Md ~ 600. Including the atomic gas component yields a total gas-to-dust 
ratio of ~ 830. Such high values are typical for IR-bright spirals and in order 
to reconcile them with the significantly lower ratio of ~ 100 obtained for the 
Milky Way a cold dust reservoir, inconspicuous at FIR wavelengths, is usually 
postulated. However, while there is good evidence for the presence of cold 
gas/dust in NGC 7469 beyond its central region, our 450 iim map and available 
interferometric ^^CO J=l-0 maps show the bright sub-mm/CO emission 
confined in the inner ~ 2.5 kpc, where a high ^^CO (J=3-2)/(J=l-0) ratio 
(~ 0.85 — 1.0) is measured. This is consistent with molecular gas at Tkm ^ 30K, 
suggesting that the bulk of the ISM in the starburst center of NGC 7469 is 
warm. Nevertheless the corresponding total gas-to-dust ratio there remains 
high, of the order of ~ 500. We argue that, rather than unaccounted cold dust 
mass, this high ratio suggests an overestimate of M(H2) from its associated 
^^CO J=l-0 hne luminosity by a factor of ~ 5 when a Milky Way value for 
this conversion is used. Finally the diffuse cold gas and dust that is the likely 
source of the observed faint extended 450 /im and ^^CO J=l-0 emission has an 
estimated total gas-to-dust ratio of ~ 50 — 160, closer to the Galactic value. 
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1. Introduction 

The role of molecular gas as the fuel of both starburst activity and an Active Galactic 
Nucleus (AGN) is now well estabhshed. In Seyfert galaxies such gas is found on scales 
ranging from the inner ~ 1 — 2 kpc "feeding" an intense circumnuclear starburst down to 
L ^100 pc from the AGN, where in the form of a geometrically thick torus obscures the 
active nucleus along certain viewing angles thus creating the difference between type 2 and 
type 1 Seyferts (Miller & Antonucci 1983; Krolik 1990 and references therein). 

Estimates of molecular gas mass in other galaxies from their velocity-integrated ^^GO 
J=l-0 luminosity Leo (in K km s^^ pc^) arc based on the so-called conversion factor 
Xco = M(H2)/Lco ~ 5 M0(K km s"^ pc^)"^ (e.g. Solomon & Barrett 1991) whose value is 
derived from studies of molecular gas in the Galactic disk. Numerous studies quantify the 
effects of the various physical conditions on this factor (e.g. Bryant & Scoville 1996; Israel 
1997). Particularly in the starburst nuclei of very luminous IR galaxies (LpiR > 10^^ Lq) 
these conditions are significantly different than in the disk environment of quiescent spirals 
like the Milky Way. More specifically in the inner ~ 1 — 2 kpc, where a starburst usually 
occurs, the gas differentiates into two distinct phases (e.g. Aalto et al. 1995) with a warm, 
diffuse and possibly non self-gravitating phase dominating the ^^GO emission. This results 
in a significant overestimate of H2 mass from Leo when a Galactic value of Xco is used 
(Solomon et al. 1997; Downes & Solomon 1998). 

The aforementioned conditions are expected in the nuclear region of the SBa Seyfert 1 
galaxy NGG 7469 (Arp 298, Mrk 1514), a luminous IR source with Lfir ^ 3 x 10^^ L© 
emanating from a powerful starburst deeply embedded into its L ~ 1 kpc circumnuclear 
region (Wilson et al. 1991; Genzel et al. 1995). High resolution ^^GO J=l-0 observations 
revealed large amounts of molecular gas within this region (Meixner et al. 1990; Tacconi 
& Genzel 1996) with mass comparable to the dynamical, and there is significant evidence 
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that application of the standard Galactic conversion factor overestimates M(H2) in the 
circumnuclear starburst of this galaxy (Genzel et al. 1995). In this paper we present 
sensitive sub-mm imaging of NGC 7469 at 850 fim and 450 fim and ^^CO J=3-2 
spectroscopy of its central region. Under the assumption of a canonical gas-to-dust ratio 
of ~ 100, our results confirm the overestimate of H2 gas mass when a Galactic value for 
Xco is used, and demonstrate the significance of CO spectroscopy and sub-mm imaging in 
offering a better assessment of the molecular gas mass and its average physical conditions. 

Throughout this work we adopt Ho = 75 km sec^^ Mpc~^ and qo = 0.5, which for 
cz=4900 km s~^ yields a luminosity distance of Dl ~ 66 Mpc for NGC 7469, where 1" 
corresponds to ~ 310 pc. 



2. Observations 

The sub-mm observations were made on the nights of 1997 December 3 and 1998 
January 15 with the Sub-mm Common User Bolometer Array (SCUBA) at the 15-m James 
Clerk Maxwell Telescope (JCMT)0. SCUBA is a dual camera system cooled to ~ 0.1 K 
allowing sensitive simultaneous observations with two arrays. The short-wavelength array 
contains 91 pixels and the long- wavelength array 37 pixels, with approximately the same 
field of view, namely ~ 2.3'. For a description of the instrument see Holland et al. (1998). 

We performed dual wavelength imaging at 450 fim and 850 /im using the 64-point 
jiggle mapping mode that allows Nyquist sampling of the field of view. ( [Holland et al 



1998|) . We employed the recommended rapid beam switching at a frequency of ~ 8 Hz and 



^The JCMT is operated by the Joint Astronomy Center in Hilo, Hawaii on behalf of 
the parent organizations PPARC in the United Kingdom, the National Research Council of 
Canada and the The Netherlands Organization for Scientific Research. 
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a beam throw of 150" in azimuth. The pointing and focus were monitored frequently using 
Uranus and CRL 618, with an expected rms pointing error of ~ 3". All maps of NGC 7469 
were "bracketed" by sky-dips that were later used to correct for atmospheric extinction. 
Conditions were generally excellent for sub-mm observations with typical opacities of 
T850 ~ 0.1 — 0.2 and T450 ~ 0.45 — 0.55 throughout our runs. 

The beam characteristics and cahbration gains (Jy beam~^) were deduced from 
an extensive archive of beam maps of Uranus and CRL 618 taken during similar periods 
as our observations and with a similar beam-throw (120"). High S/N beam maps are 
particularly important for the calibration of 450 /im images since the beam shape and the 
gain can change significantly as the dish thermally relaxes. All jiggle maps were flat-fielded, 
corrected for atmospheric extinction and edited for bad bolometers/integrations using the 
standard reduction package SURF (Jenness & Lightfoot 1998). Sky-noise was removed by 
using the bolometers in the outer ring of the two arrays (Jenness, Lightfoot & Holland 
1998) which were assumed "looking" only at sky emission. This is a good assumption for 
NGC 7469 where most of its ^^CO J=l-0 emission (and presumably most of the sub-mm 
emission from dust) lies within a ~ 4" radius (Tacconi & Genzel 1996, Tacconi et al. 2000). 

The gain and beam characteristics at 850 fim remained essentially the same throughout 
the runs. This is further corroborated by the fact that the "raw" peak and integrated 
intensities of the NGC 7469 images at this wavelength agree to within ~ 10%, which is the 
expected calibration uncertainty. Hence we co-added all the noise-weighted frames after 
subtracting a small gradient from one of them, and then scale the final one with the single 
gain factor of Ggso = 290 Jy Volt~^ beam~^. The deduced beam- width is ©hpbw ~ 

The 450 /xm beam maps revealed substantial gain variations (~ 40%) between 
observations taken early in the first shift and the ones taken later on, presumably due 
to the thermal relaxation of the dish. Thus before co-adding the images of NGC 7469 
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we scaled them by the appropriate gains of G450 = 1170 Jy Volt~^ beam~^ and 
G450 = 780 Jy Volt^^ beam^^. This scahng brought their peak and integrated intensities in 
agreement to within ~ 15%, which is comparable to the calibration uncertainty expected at 
450 /im. The average beam- width deduced for this wavelength is Ohpbw ~ 9"- 

The i^CO J=3-2 line at 345.796 GHz was observed towards the nucleus of NGC 7469 
with receiver B3 on 1999 August 6. We used the DAS spectrometer with a bandwidth 
of 920 MHz (~ 800 km s~^), and two independent channels centered on the line. The 
chop scheme employed was rapid beam switching with a 120" azimuthal throw at the 
recommended frequency of 1 Hz. The HPBW of the telescope beam at this frequency is 
~ 14" with a beam efficiency of r/mb = 0.62 (Matthews 1999). Typical system temperatures 
were of the order of Tgys ~ 600 K. After subtracting linear baselines from all spectra their 
integrated line intensities were found to agree to within ~ 10%, the typical calibration 
uncertainty for receiver B3. We then co- added all of them to produce the final spectrum. 

3. RESULTS 

The galaxy is detected as a bright sub-mm source at 850 /im where it appears 
marginally resolved, as well as at 450 /im where it clearly shows faint extended emission. 
We convolved a high-resolution map of ^^CO J=l-0 emission (Tacconi et al. 2000) from its 
original resolution of ~ 2.5" to the resolution of the 450 /im map (~ 9") and show them, 
together with the 850 /xm map, in Figure 1. 

EDITOR: PLACE FIGURE |l| HERE. 

The correspondence between the bright 450 /im and ^^CO J=l-0 emission is very good, 
suggesting that they both trace the same ISM material. The high-resolution interferometric 
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^^CO J=l-0 maps (Tacconi & Genzel 1996; Tacconi et al. 2000) show bright emission 
arising from a region with a diameter of d ~ 8", where an intense starburst is embedded 
(Wilson et ah 1991), hence the presence of warm gas and dust is expected there. Indeed 
the central region of NGC 7469, besides hosting an AGN, also harbors several supergiant 
star formation regions containing numerous (few xlO^) OB stars each (Genzel et al. 1995), 
whose intense UV light and subsequent supernova explosions will warm and disrupt the 
molecular clouds present. The warm dust of these clouds can easily dominate the global 
FIR and even the mm/sub-mm emission from this galaxy. In principle the new sub-mm 
measurements allow for a better evaluation of the dust mass and temperature. Thus in 
Table 1 we compile them together with available mm and FIR data and use them to model 
the dust emission SED, shown in Figure 2. 



EDITOR: PLACE TABLE 1 HERE 



EDITOR: PLACE FIGURE | HERE. 



The reported mm/sub-mm fluxes at 1.3 mm and 850 fim have been corrected for 
non-dust emission from CO lines, which can be rather significant (see Appendix). The fitted 
dust temperature and mass are: Td ~ 35 K and Md = 2.5 x 10^ M©. It can be seen that a 
single dust component fits the data rather well, except for wavelengths shortward of 60 fim 
where the hot dust 300 K) present in the central region (Cutri et al. 1984) makes a 
significant contribution to the total flux but contains a negligible fraction (< few x 10"'^) of 
the total dust mass. From the ^^CO J=l-0 map obtained by Tacconi et al. (2000), we find 
a total flux of Sco = (275 ± 30) Jy km s^^, which allows an estimate the global molecular 
mass from 
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M(H2) = 2.45 X 10^ Xco Sco (1) 

where D=66 Mpc and Xco ~ 5 Mq (km pc^)~^. This yields a mass of 
M(H2) ^ 1.5 X 10^° Mq and M(H2)/Md ^ 600, a typical value for IRAS galaxies 
(Young et al. 1986; Young et al. 1989; Stark et al. 1986). Including the HI gas mass 
M(HI) = 5.7 X 10^ M0 (Mirabel & Wilson 1984) yields a total gas-to-dust ratio of ~ 830. 

The presence of warm and dense molecular gas in the nucleus of NGC 7469 is 
responsible for the strong ^^CO J=3-2 line observed towards it. The high resolution ^^CO 
J=l-0 interferometer map allows for an estimate of the (J=3-2)/(J=l-0) line ratio in the 
inner d ~ 8" of this galaxy. This is possible since virtually all the bright ^^CO J=l-0 
emission lies within that region, and hence the ^^CO J=3-2 emission with its higher 
excitation requirements will have, at most, a similar size. Hence the observed main-beam 
brightness of the J=3-2 transition can be expressed as follows 

T.b=(l-e-^^) t1"^(i + Y3^c) (2) 

where t|,'"'' is the brightness temperature of the warm gas in the nucleus, x = Vln2 d/GnpEW 
(a disk source assumed), and C = T^.'^^T^'^^ is the ratio of radiation temperatures for J=3-2 
of the cold and warm gas phase. The former dominates at larger galactocentric distances in 
this galaxy (Papadopoulos & Seaquist 1998), and fills the rest of the 14" beam. 

Substituting d = 8" and Ghpbw = 14" we obtain t{7^ = 4.94 (1 + 3.94C)"^ T„,b, which 
we plot in Figure 3 (for C = 0) together with the ^^CO J=l-0 brightness temperature 
averaged over the inner diameter of ~ 8". 



EDITOR: PLACE FIGURE | HERE. 
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The excellent agreement between the two spectral line profiles demonstrates that 
they arise from the same region, which the high-resolution -"^^CO J=l-0 maps (e.g. 
Tacconi & Genzel 1996) allow us to identify with the inner 8". The velocity-averaged 
1^32 = (3 — 2)/(l — 0) ratio for that region is R32 ~ 1.1 (C=0). If the typical physical 
characteristics of the cold phase, namely n(H2) 3 x 10^ cm~^ and Tkin fa 10 K 
(Papadopoulos & Seaquist 1998), dominate beyond the central 8" while warm and dense gas 
(n(H2) > 10^ cm"^, Tkm > 20 K) is present in the central region, then C ~ 0.030 - 0.045. 
This yields a range R32 ~ 0.95 — 1.0, which for a gaussian source brightness will be 
somewhat lower, namely ~ 0.80 — 0.85. 

Such high R32 ratios are characteristic of an optically thick and thermahsed J=3-2 
transition at gas temperatures of Tkin ^ 30 K, consistent with the deduced dust temperature 
of ~ 35 K. A simple Large Velocity Gradient (LVG) analysis (e.g. Richardson 1985) yields 
a variety of conditions that can produce R32 ~ 0.8 — 1.1, all of them characterized by 
Tkin <; 30 K and with the most typical temperature being ~ 50 K. This suggests that most 
of the molecular ISM in the starburst nucleus of NGC 7469 is indeed warm with its dust 
emission dominating the global mm/sub-mm/FIR spectrum (Figure 2). 



4. Discussion 

The high R32 ratio and gas temperature inferred for the central region of NGC 7469 are 
in sharp contrast with the low R21 = (2 — 1)/(1 — 0) ~ 0.5 and Tkin ~ 10 K characterizing 
the global CO emission from its disk (Papadopoulos & Seaquist 1998). Such steep excitation 
gradients in spiral galaxies, particularly the ones with starburst nuclei, are expected and 
known to exist in numerous cases (Knapp et al 1980; Wall & Jaffe 1990; Wall et al. 1991; 
Eckart et al 1991; Harris et al. 1991; Wild et al. 1992; Aalto et al. 1995). This is further 
corroborated by the new generation of mm and sub-mm bolometers that allow sensitive 
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imaging of the dust emission out to large galactocentric distances and show large masses of 
cold dust (Td = 10 — 15 K) residing in spiral disks away from the nucleus (e.g. Neininger et 
al. 1996; Dumke et al. 1997; Alton et al. 1998; Papadopoulos & Seaquist 1999a). 

It is often argued that such a cold dust component is responsible for a systematic 
underestimate of dust mass when only FIR data are used. This would naturally lead to the 
large M(H2)/Md ~ 500 - 600 ratios (M(H2 + HI) ~ 1000 when HI is included) found for 
spiral galaxies in contrast to ~ 100 — 150 found for the Milky Way (e.g. Devereux & Young 
1990 and references therein). On the other hand extragalactic M(H2) estimates are thought 
to be accurate to within a factor of ~ 2 (Young & Scoville 1991), and particularly for spirals 
to within ±30% (Devereux & Young 1990). However recent sensitive sub- mm imaging at 
450 /xm and 850 fj,m of IR-bright spirals hke NGC 891 (Alton et al. 1998) and NGC 1068 
(Papadopoulos & Seaquist 1999a) reveals that the warm and cold dust are spatially well 
separated, with the former being concentrated mainly in the inner ^ 2 kpc where significant 
star formation occurs and the latter residing at larger galactocentric radii. 

The case of NGC 1068 is particularly relevant since its host galaxy properties are 
remarkably similar to NGC 7469 (Wilson et al. 1991), with both galaxies harboring a 
central starburst and having similar FIR, ^^CO J=l-0 and HI luminosities. The relative 
proximity of NGC 1068 allowed extensive multi-line CO and sub-mm imaging of its inner 
2.5 kpc (Papadopoulos & Seaquist 1999a, 1999b) which showed that the bulk of the gas and 
dust in that region is warm with M(H2)/Md ~ 330. The latter is essentially identical to the 
total gas-to-dust ratio since most of the gas in the circumnuclear starburst of NGC 1068 is 
molecular. This ratio is actually ~ 30 — 40% higher still if the bright ^^CO J=3-2 emission 
is subtracted from the 850 //m flux used to estimate the dust mass. In the regions beyond 
the starburst nucleus of NGC 1068 emission from cold dust (Td ~ 10 — 15 K) dominates the 
sub-mm bands, while low- brightness ^^CO J=l-0 and bright HI emission reveal the diffuse 
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H2 and HI that make up the gas reservoir (Papadopoulos & Seaquist 1999a). Interestingly 
the total gas-to-dust ratio in those regions is ~ 70 — 150, close to the Galactic value. 

The larger distance of NGC 7469 precludes a similar detailed study, however the 
high-resolution ^^CO J=l-0 and 450 /xm maps (Figure 1) together with the lower limit 
of ~ 35 K for the dust/gas temperature in its central 8" allow for a firm lower limit 
on the gas-to-dust ratio in that region. We measure S45o(r < 4") = (0.80 ± 0.08) Jy, 
probably a slight overestimate owing to the larger beam area at this wavelength coupling 
to a region somewhat larger than the inner 8", beyond which only faint sub-mm 
emission from cold dust is expected. The velocity-integrated ^^CO J=l-0 line flux is 
Sco(r < 4") = (175 ± 25) Jy km s~^, estimated from the original high resolution map. 
Assuming that H2 dominates the gas phase in the starburst region, a Galactic Xco, and 
Td = 35 K, yields a total gas-to-dust ratio of ~ 500. 

Hence for both NGC 7469 and NGC 1068 it can be convincingly argued that the high 
gas-to-dust ratios found for their IR-luminous starbursting central regions are not due to 
the presence of significant amounts of cold dust but to an overestimate of the molecular gas 
mass in such environments. In the case of NGC 7469 this conclusion is further supported by 
a dynamical study of its central region concluding that Xco is ~ 1/5 of the Galactic value 
(Genzel et al. 1995). This would bring the gas-to-dust ratio estimated for that region in 
good accord with the Milky Way value with no need for an unaccounted mass of cold dust. 

4.1. The systematic overestimate of M(H2) in stairburst environments 

There has been mounting evidence that the Galactic value of Xco overestimates 
the molecular gas mass in the intense starburst environments of luminous IR galaxies 
(Lfir > 10^^ Lq by a factor of ~ 5 (Solomon et al. 1997; Downes & Solomon 1998). This 
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seems to be due to a two-phase differentiation that the molecular gas undergoes in such 
environments (Aalto et al. 1995; Downes & Solomon 1998). The phase that dominates the 
^^CO emission is diffuse, warm and possibly non self-gravitating, the latter being the main 
reason for the overestimate of molecular gas when a Galactic value for Xco is used. 

Indeed a standard expression for this conversion factor is (e.g. Bryant & Scoville 1996) 

Xco = 2.1 ^ Q Mo (K km s-')-' (3) 

where n (cm~^) is the average gas density and Tj, (K) is the average brightness temperature 
of ^^CO J=l-0 of the molecular cloud ensemble. The factor Q = SYyir/SY accounts for the 
non-virial hnewidth of the "average" cloud, for self gravitating clouds it is Q = 1. 

For typical conditions in the Galaxy, namely Q fa 1, n 300 cm~^ and Tkin ~ 15 K 
(Tb 8 K), one obtains Xco ~ 5 Mq (K km s~^)~^, the typical Galactic value. In the 
diffuse molecular phase that is present in starburst nuclei, it can be Q < 1 (e.g. Solomon et 
al. 1997), thus a Galactic Xco will overestimate the molecular gas mass present. Moreover, 
even in the case of self-gravitating clouds, a diffuse (n 10^ cm~^) but warm phase 
(Tkin ~ 60 K) dominating the ^^CO J=l-0 emission will have Tb = 30 - 40 K (no ~ 2 - 4) 
and thus yield a Xco factor that is ~ 2 — 3 times smaller than the Galactic value. 

If such a gas phase is characteristic of the starburst regions in luminous IR galaxies, 
then a systematic overestimate of H2 mass rather than an underestimate of dust mass is 
responsible for the high M(H2)/Md ratios found for them. Because the effect is systematic 
it can produce both a high total gas-to-dust ratio and a small dispersion around its mean. 
Hence earlier claims that such small dispersion (~ ±30%) observed in IR-luminous spirals 
(Devereux & Young 1990) argues in favor of a similar accuracy in molecular gas mass 
estimates are not supported by this picture. 
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4.2. The cold gas and dust in NGC 7469 

The global excitation of CO in NGC 7469 seems to be dominated by cold and 
diffuse molecular gas out to a scale of d ^ 20" (Papadopoulos & Seaquist 1998). A 
gaussian fit of the 850 /xm emission (Figure 1) yields a source size of ~ 17.2" x 18.5" 
(FWHM), clearly larger than the beam at this wavelength. Considering the geometric 
mean of ©s ~ 18" to be the observed diameter of the source (assumed to be a disk), 
and for a gaussian beam with ©hpbw = 15", we obtain an intrinsic source diameter of 
d = [2 (ln2)-i (©2 - ©Ipbw)]^^^ ~ 17"- This is ~ 2 times larger than the size of the bright 
^^CO J=l-0 emission and comparable to the source size seen at 450 //m (Figure 1). 

In NGC 1068 the warm gas/dust resides in the inner d ~ (2.7 — 3.4) kpc, which 
is comparable to the size of a similar region in NGC 7469, both regions containing an 
embedded starburst. Furthermore, given the similarities of the disk properties between 
these two galaxies, it is reasonable to identify the extended faint sub-mm, ^^CO J=l-0 
emission in NGC 7469 with emission from cold dust/gas. Indeed simply scahng the total 
size of the faint sub-mm emission from cold dust in NGC 1068 (~ 6 kpc) to the distance of 
NGC 7469 yields an angular size of ~ 20", closely matching the one observed for the latter. 

It is intriguing that faint ^^CO J=l-0 emission is detected out to significantly larger 
radii than the inner radius of r = 4" (Figure 1). Assuming that, as in NGC 1068, this 
emission "marks" cold and diffuse H2 coexisting with cold dust (Ta = 10 — 15 K) and the 
bulk of HI, yields a ratio [M(H2 + HI)]/Md ~ 50 - 160, close to the Milky Way value. 

4.3. Warm versus cold gas/dust: sub-mm imaging and CO spectroscopy 

It is striking that in the spatially integrated dust emission of this galaxy the cold 
dust is inconspicuous even after the inclusion of the sub-mm data. On the other hand its 
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molecular gas counterpart dominates the global excitation characteristics as revealed by the 
low ^^CO (J=2-l)/(J=l-0) line ratio (Papadopoulos & Seaquist 1998). 

This could simply be due to the fact that spectral line luminosity is sensitive to the 
local gas density as well as its temperature and total mass but only the latter two determine 
the continuum luminosity from dust (for a fixed gas-to-dust ratio). In other words two 
molecular clouds with the same mass and gas/dust temperature will emit the same dust 
continuum (optically thin case) but their CO line fluxes can differ drastically if their average 
H2 densities are different. This will occur mainly in the density regime of sub-thermal 
excitation of the observed CO line(s). 

In NGC 7469 the relative amounts of warm and cold dust can be estimated from the 
relation 



m = — j-^ 
Mi") 



g32/Tc _ I 
e32/Tw _ 1 



c(c) 
£450 

q(w) ' 



^450 

where (w) and (c) denote the quantities corresponding to the warm and cold dust 
respectively. It is s2 = S45o(r < 4") = 0.8 Jy and Sgo = 8450(4" < r < 11") = 0.5 Jy. 
Hence for T^ ;^ 35 K and Tc = 10 - 15 K we obtain m ;^ 3 - 10. 

Even for m ~ 10, the spatially averaged sub-mm emission will not necessarily 
reveal the presence of the cold component. Indeed, since the ratio r = S450/S850 is 
temperature-sensitive when Td ^ 30 K, it is expected to differ if a cold dust component is 
present along with the warm one. Assuming spatial averaging of the emission from the two 
dust components the observed r can be written as follows 



1 + m f45o(T 
1 -Km f85o(T^,Tc)_ 

where r^ is the ratio corresponding to warm dust alone and. 



(5) 
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gT;\/Tw ^ ll C 

For T„ = 40 K, Tc = 10 K and m = 10 we obtain r f« 0.70 x r^, which is barely 
discernible from r^ given that the value of r carries a ~ 20% uncertainty due to the 
cahbration uncertainties of SCUBA at 450 iim and 850 iim. 

On the other hand besides temperature the spectral line ratios are also sensitive to the 
local gas density. Thus a cold and sub-thermally excited molecular gas phase may dominate 
the observed global line ratio even in the presence of a warm and thermalized phase. For 

example, the global R21 = (2 — 1)/(1 — 0) ratio can be expressed as 

l + mx(RgMT^) l ... 
1+mx J' ^' 

where x = X^q/X^q is the ratio of the CO(l-0)-to-H2 conversion factors for the two gas 
phases. Assuming x ~ 1/5, m ~ 10 and, R^^ ~ 1, we obtain R21 ~ 1/3 (1 + 2 R^^^). 

In NGC 7469 it is R21 ~ 0.5 (Papadopoulos & Seaquist 1998), hence the cold- phase 
ratio will be R^i 0.25, and for Tkin = 10 K, the latter corresponds to densities of 
n(H2) ~ 10^ cm~^. For the same temperature but thermalized and optically thick ^^CO 
J=2-l it is R21 ~ 0.8, yielding R21 ~ 0.85 which is significantly larger than observed. 

The aforementioned simple analysis demonstrates the significance of combining sub-mm 
and CO observations in evaluating the physical conditions and the mass of molecular gas. 
The sub-mm measurements offer an independent means of estimating gas mass under the 
assumption of a canonical gas-to-dust ratio of ~ 100. For H2-dominated gas this mass 
can then be compared to the one deduced from the ^^CO J=l-0 luminosity and the Xco 
conversion factor and hence the influence of the physical conditions on Xco can be assessed. 
On the other hand, if a cold and diffuse molecular gas phase is present along with a warm 
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one but only global averages are available, CO line ratios can be more sensitive to the 
presence of the cold ISM phase than sub-mm intensity ratios. 

5. Conclusions 

We presented sensitive 850 fim, 450 fim maps of the Scyfert 1 galaxy NGC 7469 and 
a measurement of the ^^CO J— 3-2 line towards its starbursting central region. Our main 
conclusions can be summarized as follows 

1. The FIR/sub- mm/mm spectrum of this source is dominated by a warm dust 
component with a temperature of ~ 35 K and a global molecular gas-to-dust ratio of 
~ 600, both typical for IR-luminous spirals. Including the atomic gas mass yields a total 
gas-to-dust ratio of ~ 800. 

2. The warm dust and gas hes in the inner 8" (~2.5 kpc) where a star burst is 
embedded. In this region we find no evidence for a significant mass of cold dust, yet the 
ratio of the mainly molecular gas to the dust mass is ~ 500, still five times larger than 
the Galactic value. We argue that this is the result of a systematic overestimate of H2 
mass by a factor of ~ 5 when a Galactic value for Xco = M(H2)/Lco is used in starburst 
environments. 

3. On larger scales (radius ^ 1.2 kpc) the ISM in NGC 7469 is dominated by cold, 
sub-thermally excited gas where the faint 450 /im and ^^CO J=l-0 emission originate. The 
gas-to-dust ratio for this phase, is ~ 50 — 150; in better accord with the Milky Way value. 
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A. Spectral line contribution to mm/sub-mm bands 



In mm and sub-mm continuum observations one must correct for non-dust emission 
contributions. The most significant comes from molecular spectral lines and can reach up 
to ~ 60% of the total flux observed with a typical bolometer bandwidth (Gordon 1995). 
In warm cores of Orion such line contributions are found to be of the order of ~ 30% in 
850 //m and 450 //m bands (Johnstone & Bally 1999 and references therein). However 
in the extragalactic domain a typical beam of a mm/sub-mm telescope samples areas of 
several hundred parsecs, and the contributions from warm cores (< 1 pc) in star forming 
regions are negligible because of the resulting spatial dilution. Over such scales only the 
emission from the three lowest rotational transitions of ^^CO can be ubiquitous and bright, 
since these transitions can be easily excited in the general conditions prevailing in the ISM 
in quiescent or starburst environments. 

In the case of NGC 7469 we expect that all the signiflcant spectral line emission from 
the ^2 CO J=2-l (1.3 mm band) and ^^CO J=3-2 (850 //m band) lines arises from the 
region of bright ^^CO J=l-0 emission in the central 8" where the starburst resides. In the 
case of the 1.3 mm band only line- free channels were used to produce the continuum map 
(Tacconi, private communication). The 850 /im band includes the ^^CO J— 3-2 line whose 
contribution can be estimated from our observation of this line towards the nucleus, namely 



where Imb (K km s~^) is the velocity-integrated main-beam brightness of ^^CO J=3-2, 
I'o = 345.796 GHz, and fimb is the gaussian beam area at this frequency. The term in the 
parenthesis corrects for the beam-source geometrical coupling assuming a disk source of 
diameter d, where x = -\/ln2 d/©HPBW, and Aub is the bolometer bandwidth. 




(Al) 
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For i^CO J=3— 2 in NGC 7469 we measured I^b = (77 ±8) K km s'^ (Figure 3), hence 
for Ghpbw = 14", Az/b = 30 GHz (Holland, private communication) and d = 8" we obtain 
a line contribution of AS = 63 mJy, which amounts to ~ 40% of the total flux at 850 yum. 
The ^^CO J=3-2 is expected to have ^ 0.1 of the ^^CO J=3-2 flux, and hence contribute 
^ 4% of the total flux in 850 /xm. 

In the case of the 450 /im band the only CO hue that may contribute to the observed 
flux is ^^CO J=6-5 at 661.067 GHz. However we find that this hne, with its high excitation 
density (n > 10^ cm^^), will not contribute significantly to the 450 //.m band (i.e ^ 2%) 
over the spatial scales relevant here, even in a starburst environment. 
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Fig. 1. — Top: the 850 /im map with a FWHM= 15" beam shown at the lower left, the 
contours are (—3, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30) x agso, with o-gso = 5 mJy beam~^. 
Bottom: the 450 /xm map (contours) overlaid to the integrated ^^CO J=l-0 emission (grey 
scale) at a common resolution of 9" (FWHM beam shown at the bottom left). The contours 
are (—4, 4, 6, 8, 10, 14, 16, 18) x (T450, with (J450 = 40 mJy beam~^ and the grey scale is 1=6- 
200 Jy beam~^ km s^^ (o'(I) = 1 Jy beam~^ km s~^) 
The map center is at RA: 23^^ 03°^ 15^6, Dec: +08° 52' 26" (J2000). 

Fig. 2.— A x^-fitted SED for the mm/sub-mm/FIR data of NGC 7469 (Table 1). An 
optically thin isothermal dust reservoir has been assumed with an emissivity law of a = 2 
and emissivity at 200 //m (1196 GHz) of ko = 10 cm^ gr"^ (Hildebrand 1983). The flux at 
25 iim is not used in the fit. 

Fig. 3.— The ^^CO J=l-0 and J=3-2 brightness temperatures for the inner 8" of NGC 7469 
(see text). The thermal rms uncertainty across the band is STj-j^si^ — 0) = 30 mK, and 
S%^s{^ - 2) = 40 mK. 
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Table 1 

NGC 7469: The mm, sub-mm and FIR data 



Wavelength 


Flux 


Reference 


1.3 mm 


19 ± 5 mJy 


Tacconi et al. 1999 " 


850 urn 


91 ± 13 mJy 


this work ^ 


450 yum 


1.30 ±0.14 Jy 


this work ^ 


100 urn 


34.9 ± 0.6 Jy 


Soifer et al. 1989 


60 /im 


27.7 ± 0.04 Jy 


u u 


25 /im 


5.84 ± 0.05 Jy 


u u 



^ Derived from the 1.3 mm continuum map for a radius of r = 4" that contains all the mm 
emission and corrected for non-thermal emission using the data from Wilson et al. 1991. 

^ Fluxes are estimated for a radius of r = 11" and the appropriate error-beam corrections 
(Sandell 1997) are fgso = 0.80, f45o = 1.05 — 1.10, derived from Uranus and CRL 618 maps. 
The 850 iim flux has also been corrected for contribution from the ^^CO J=3-2 line (see 
Appendix). The errors reported include (besides the thermal rms errors) a ~10% (850 //m) 
and a ~15% (450 /xm) calibration uncertainty. 
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